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We present optically tunable magnetic 3D metamaterials at terahertz (THz) frequencies which 
exhibit a tuning range of ^30% of the resonance frequency. This is accomplished by fabricating 3D 
array structures consisting of double-split-ring resonators (DSRRs) on silicon-on-sapphire, fabricated 
using multilayer electroplating. Photoexcitation of free carriers in the silicon within the capacitive 
region of the DSRR results in a red-shift of the resonant frequency from 1.74 THz to 1.16 THz. 
The observed frequency shift leads to a transition from a magnetic-to-bianisotropic response as ver- 
ified through electromagnetic simulations and parameter retrieval. Our approach extends dynamic 
metamaterial tuning to magnetic control, and may find applications in switching and modulation, 
polarization control, or tunable perfect absorbers. 



The advent of artificial electromagnetic materials has 
provided unique routes to engineer photonic media, lead- 
ing to the realization of intriguing phenomena. This in- 
cludes the successful demonstration of invisibility cloak- 
ing [T-'S and negative refraction [^J |5j . Scale invariance 
of the underlying equations enables translation of phe- 
nomena realized in one region of the electromagnetic 
spectrum to others. Metamaterials and in particular, 
split-ring resonators (SRR), can couple to the electric, 
magnetic or both fields of incident electromagnetic radi- 
ation. This leads to a Lorentzian-like effective response 
in the permittivity, permeability, and bianisotropy, de- 
scribed by the general constitutive relations as: 

/ p \ ^ / aEE C^EH \ f \ 
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where aEE {<^hh) is the tensorial purely electric (mag- 
netic) polarizability, quantifying the electric (magnetic) 
dipole induced purely by a incident electric (magnetic) 
field and the off-diagonal aEH {(^he) is the tensorial 
magnetoelectric coupled polarizability, quantifying the 
electric (magnetic) dipole induced by an incident mag- 
netic (electric) field. The Onsager relations require 
aEH = ~o^JJ^ [6^. The bianisotropic parameter ^ relates 
to the magnetoelectric polarizability by ^ = icoaEH^ 
where cq is the velocity of light in free space. 

By tailoring the geometry or configuration of subwave- 
length metallic unit cells, the effective parameters can be 
specified and controlled, leading to myriad routes to con- 
trol the amplitude, frequency, and phase of the incident 
electromagnetic radiation. For example, by indepen- 
dently tuning the electric and magnetic resonance of mul- 
tilayer composites, it is possible to match the impedance 
to free space forming a perfect absorber with minimized 
reflectance [7]-[9]. This ability to manipulate incident ra- 
diation is particularly important at terahertz frequencies 
given the relative paucity of devices. Recent progress on 
tunable metamaterial through voltage [TOl [TT], optical 
[T2UT4] . thermal [T5HT7] . and mechanical 58l[19] tuning. 
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FIG. 1: (Color online) Schematic of the design of the broad- 
band tunable three dimensional metamaterials at THz fre- 
quencies and scanning electron micrographs (SEM) of fabri- 
cated 3D structures, (a) Schematic principle of the photoex- 
citation tuning of the 3D metamaterials. (b) Typical dimen- 
sions for unit cell of DSRRs: p = 50 /xm, h = 33 /xm, g = 5/im, 
and 1 = 36 /xm. (c) The oblique close-up view of fabricated 
unit cell of DSRRs in copper standing on a sapphire substrate, 
(d) Oblique view of arrayed hybrid 3D metamaterials. 



has shown its potential to fill the "THz gap" . 

The majority of previous studies are based on planar 
metamaterials, and therefore result in tuning the electric 
response. It is desirable to explore new configurations to 
access magnetic control of metamaterials which will, in 
turn, provide additional flexibility in tailoring the electro- 
magnetic response for potential applications such as tun- 
able absorbers or polarizers. Recently, three dimensional 
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FIG. 2: (Color online) (a) Experimentally measured fre- 
quency dependent THz electric field transmission. The re- 
sponse is shown for increasing laser excitation power, (b) Full 
wave electromagnetic simulation results of the transmission 
versus frequency for various values of the silicon conductivity. 
The dashed black curve is a simulation with the bottom gap 
shorted with copper. 



metamaterials at terahertz frequencies have been demon- 
strated with the successful demonstration of negative re- 
fraction [20 and negative permeability [21 , with one ex- 
ample of a tunable three dimensional (3D) metamate- 
rial [22] . We demonstrate optically tunable magnetic 3D 
metamaterials exhibiting a tuning range of ~30% of the 
resonance frequency. This is accomplished by fabricating 
3D array structures consisting of double-split-ring res- 
onators (DSRRs) on silicon-on-sapphire, fabricated using 
multilayer electroplating. Photoexcitation of free carriers 
in the silicon within the capacitive region of the DSRR 
results in a red-shift of the resonant frequency from 1.74 
THz to 1.16 THz. Further, the observed frequency shift 
results in a transition from a magnetic-to-bianisotropic 
response as verified through electromagnetic simulations 
and parameter retrieval. 

The tuning principle and unit cell of the three dimen- 
sional metamaterial structure are depicted in Fig. [l] (a) 
and (b). The arrayed 3D structures consist of stand- 
up DSRRs with gaps on the top and bottom, forming 
two capacitors connected in series in an equivalent LC 
circuit. Silicon, the photoactive material in these struc- 



tures, is incorporated into the bottom gap of the DSRR. 
This structure is resonant when the THz wave is normally 
incident on the sample with the E field parallel to the 
sides with the gaps (H then passes through the DSRRs). 
In the absence of photoexcitation, the THz pulses cou- 
ple to the resonators magnetically, inducing a circulat- 
ing current in the ring. In particular, the symmetry is 
such that the electric field does not couple to the DSRRs 
(except for substrate-induced symmetry breaking as dis- 
cussed below). Upon optical excitation of the silicon, free 
carriers are generated resulting in an increased conduc- 
tivity. This decreases the capacitance of the lower gap 
which, in turn, increases the total capacitance resulting 
in a redshift of the LC resonance. With sufficient pump 
power, the bottom capacitor can be effectively shorted 
leading to the maximum possible frequency shift. Im- 
portantly, with the lower gap shorted, the symmetry of 
the DSRR is reduced such that the incident E field now 
drives circulating currents in the resonator (in addition 
to the H field). That is, photoexcitation induces intrinsic 
bianisotropy. 

Multilayer electroplating (MLEP) in combination with 
conventional optical lithography is utilized to fabri- 
cate the 3D metamaterials on a commercial silicon-on- 
sapphire (SOS) wafer consisting of 600-nm (100) intrin- 
sic silicon on top of a 530 /im-thick R-plane sapphire 
substrate. The fabrication of 3D DSRRs begins with 
patterning of the distributed silicon pads using reactive 
ion etching (RIE). Then two 1.5-/im-thick bottom cop- 
per laminas were electroplated on a copper seed layer, 
forming a 5-/im-wide gap between the two laminas. The 
laminas were patterned along the primary flat orientation 
of the wafer, which is 45^ counterclockwise from the pro- 
jection of the C-axis on the R-plane. Thus, the refractive 
indices on two axes are the same. The patterned silicon 
pads sit beneath the gap. After removing the photore- 
sist, a 30- /im thick layer of photoresist (PR) AZ9260 was 
patterned with two 6-/im diameter holes at both ends of 
the laminas to enable electroplating. Without removing 
the thick PR, another 100-nm seed layer of copper was 
evaporated on the top without any wrinkling. After one 
thin layer PR was patterned, 1.5- jam laminas that are 
the same as the bottom ones were fabricated to connect 
the standing pillars. As a final releasing step, the PR 
was removed with acetone, and the copper seed layers 
are etched away in dilute acid. Using MLEP, the height 
of the DSSRs can be determined by the thickness of PR, 
which is controlled by the spinning rate. Fig. [TJc) and 
(d) show the SEM images of the fabricated 3D metama- 
terials sample. 

THz-TDS was used to characterize the frequency- 
dependent electromagnetic response of the photoactive 
3D-DSRRs. A bare sapphire substrate with the same 
thickness as the sample substrate was used as the ref- 
erence. The sample and reference were mounted on the 
holder with the same orientation as to minimize sapphire 
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birefringence. The THz pulses were at normal incidence 
to the substrate with E-field parallel to the side of the 
DSRRs with the gaps (see Fig. 1(a)). Photoexcitation 
of free carriers in the silicon was achieved with 1.55-eV, 
35-fs ultrafast pulses at a repetition rate of 1 kHz. The 
optical pump pulse was set to arrive 10 ps before the THz 
probe beam to achieve a near steady-state accumulation 
of carriers due to their long lifetime in silicon. To ensure 
homogeneous excitation, the pump laser beam was 8 mm 
in diameter while the incident THz had a 3 mm diameter. 

The experimental results are shown in Fig. |2|a). The 
black curve shows the response without photoexcitation. 
There is a very sharp resonance at 1.74 THz (cji) with 
a minimum transmission of 35%. As the pump power is 
increased to 450 mW, a pronounced resonance shift to 

I. 16 THz {UJ2) is observed with a transmission of 60%. 
Further, a tunability of 30% of the resonance frequency 
is achieved. At higher powers, a large density of carriers 
result in significant carrier-carrier scattering, saturating 
the conductivity. For example, increasing the excitation 
power from 450 mW to 900 mW) corresponding to a 
fiuence of 1.8 mJ cm~^) leads to only a minor change 
in the DSRR electromagnetic response. 

Full wave electromagnetic simulations were performed 
using CST microwave studio. To compare with experi- 
ments, the conductivity of the silicon comprising the sub- 
strate capacitor was varied in the simulations. The pho- 
toconductive silicon is simulated as dielectric with Cgi = 

II. 7 and a pump-power dependent conductivity of a si- 
The substrate sapphire was modeled as lossless dielec- 
tric with Csapphire = 10.5. The copper was modeled as 
a lossy metal with a frequency independent conductiv- 
ity of 5x10^ {Qcm)~^^ based on four-point probe mea- 
surements of an electroplated continuous Cu film. The 
solid black curve in Fig. 3(b) corresponds to mW 
pump power with asi taken as 0.01 {Qcm)~^, while for 
illumination of 900 mW, the conductivity is set as 600 
{Qcm)~^. As shown in Fig. [2|b), the simulated and 
measured frequencies match each other quite well, high- 
lighting that the observed changes in the electromagnetic 
response arise from changes in the Si conductivity. 

For our experiments, the electric field is normal to the 
gaps and the magnetic field is normal to the plane of 
rings. As such, the electric and magnetic components of 
the THz waves can, in principle, couple to the DSRRs 
leading to a bianisotropic response. However, without 
photoexcitation (mode 1), inversion symmetry of the ring 
(in the propagation direction) eliminates the possibility 
of electric field excitation resulting in a pure mangetic re- 
sponse (the substrate breaks this inversion symmetry and 
leads to measurable, but small, bianisotropy - see below 
and Fig. 3). Following photoexcitation, the shorted bot- 
tom gap breaks the inversion symmetry, leading to strong 
electric field excitation and a bianistropic response. 

To further analyze the electromagnetic response of 
modes 1 and 2, the real and imaginary parts of e, /i and 
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FIG. 3: (Color online) Retrieved constitutive parameters of 
(a) permeability /x, (b) bianisotropy parameter ^, and (c) per- 
mittivity e, from full wave electromagnetic simulations of the 
3D DSRRs. The blue lines represent the parameters at mode 
1 (i.e. no photoexcitation). The red lines represent the pa- 
rameters at mode 2 assuming the bottom gap is shorted by 
replacing the silicon by copper. Solid lines are real parts and 
dashed lines are imaginary parts. 



^ were extracted from simulations using a bianisotropic 
retrieval method [23H25] . As shown in Fig. |3jb), at mode 
1, the permeability exhibits a strong Lorentzian-like re- 
sponse at the resonance yielding negative values in a cer- 
tain span. However, as mentioned above, the substrate 
breaks the symmetry and induces a small bianisotropy as 
the blue curves in Fig. [s] (b) indicate Hence, there 
is also a small resonant response for the permittivity at 
mode 1. When the bottom gap is totally shorted (mode 
2), the 3D DSRR array becomes fully bianisotropic with 
both the electric and magnetic field resonantly driving 
the array. The bianisotropic response is now dominated 
by the DSRRs. Although the permittivity and perme- 
ability present negative values, the refractive index is 
still positive due to the strong bianisotropy. Nonethe- 
less, our metamaterials can be used to effectively switch 
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FIG. 4: (Color online) Simulation of the distribution of sur- 
face current density and electric field in the direction of nor- 
mal to the gaps at two modes. Left: Surface current and 
electric field at mode 1 with a magnetic dipole of m and elec- 
tric dipoles of pi and p2. Right: Surface current and electric 
field at mode 1 with a magnetic dipole of vn and an electric 
dipoles of pi . The electric field and the surface current are 90 
degrees out of phase. The colorbars represent the simulated 
excited field strength relative to the incident field. 



In summary, we have successfully fabricated and char- 
acterized broadband tunable 3D hybrid metamateirals at 
THz frequencies. With photoexcitation of silicon, we ob- 
served over 30% redshift of the resonant frequency, in 
agreement with numerical simulations. Parameters re- 
trieval indicates that our 3D structure can be effectively 
switched from a mode with a strong magnetic response to 
a fully bianisotropic response. Our tunable 3D metama- 
terials can be potentially used as a tunable notch filter, 
high speed modulator. With proper design on the phase 
shift between two modes, artificial polarizers can also be 
realized. 

We acknowledge partial support from DOD/Army Re- 
search Laboratory under Contract No. W911NF-06-2- 
0040, NSF under Contract No. ECCS 0802036, and 
AFOSR under Contract No. FA9550-09-1- 0708. The au- 
thors would like to thank the Photonics Center at Boston 
University for all the technical support throughout the 
course of this research. 



the DSRRs with the electromagnetic response changing 
from near purely magnetic to bianisotropic. 

From the measurements and simulations, the mode 2 
resonance is broadened in comparison to mode 1. While 
there are losses from excited carriers in the silicon, this 
is not the dominant source of the broadening. This is 
evident in Fig. 2(b) showing that the simulation with 
the lower gap shorted with copper (black dashed line) 
is still considerably broader than the mode 1 resonance. 
In fact, the broadening of mode 2 results from increased 
radiation damping [27 . Additional insight into the elec- 
tromagnetic response and the nature of this damping is 
obtained through consideration of the local electric field 
and current density (from simulations) at the unit cell 
level. Fig. [4] shows the simulated surface currents and 
the electric field normal to the gap at both modes re- 
spectively. We note that the peak of the surface cur- 
rents is 90 degree out of phase with the peak of electric 
field. For mode 1, the circulating current causes charge 
accumulation at both gaps forming two electric dipoles 
Pi and p2 in opposite directions, which can be equiva- 
lently taken as a quadrupole as shown in the left side 
of Fig. 4. As such, radiation damping of the DSRRs 
for mode 1 is at the level of magnetic dipole and electric 
quadrupole terms in a multipole expansion of the elec- 
trodynamic response. In contrast, for mode 2 (in which 
the bottom gap is effectively shorted), there is only the 
top gap, which is equivalent to an electric dipole pi (right 
side. Fig. 4). Consequently, radiation damping for mode 
2 is dominated by the electric dipole contribution, lead- 
ing to a decrease in the quality factor of the resonance. It 
is important to note that this is intrinsic to the SRRs and 
is not a spurious contribution that can, for our geometry, 
be minimized. 
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